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Abstract 



> 
O 

The theory of weak measurement, proposed by Aharonov and coworkers, has 
been applied by Steinberg to the long-discussed traversal time problem. The 
uncertainty and ambiguity that characterize this concept from the perspective 
of von Neumann measurement theory apparently vanish, and joint probabil- 
ities and conditional averages become meaningful concepts. We express the 
Larmor clock and some other well-known methods in the weak measurement 
formalism. We also propose a method to determine higher moments of the 
traversal time distribution in terms of the outcome of a gedanken experi- 
ment, by introducing an appropriate operator. Since the weak measurement 
approach can sometimes lead to unphysical results, for example average neg- 
■ ative reflection times and higher moments, the interpretation of the results 

^h, obtained remains an open problem. 

3 ' I. INTRODUCTION 

cr 

In the last few years a new approach to measurement in quantum mechanics has been 
developed by Aharonov and coworkers |],|2[]. Their "weak measurement" approach differs 
from the standard one (formalized by von Neumann |3[]) in that the interaction between the 
measuring apparatus and the measured system is too weak to trigger a collapse of the wave 
function. Although an individual weak measurement of an observable has no meaning, one 
can obtain the expectation value to any desired accuracy by averaging a sufficiently large 
number of such individual results. 

Avoiding wave function collapse allows the simultaneous measurement of non-commuting 
observables (no violation of the uncertainty principle occurs because the individual measure- 
ments of each observable are very imprecise). It also allows a sound definition of conditional 
probabilities and their distribution: since the system evolves after the measurement as if 
unperturbed, it is possible to define averages of a quantity conditioned to a given final state 
of the system. Moreover - and this point is important if we are interested in the duration 
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of some process - a typical weak measurement is extended in time, i.e., the interaction be- 
tween the meter and the system is not impulsive, but has a finite duration. As Steinberg 
has shown, [fflpfl all these features make weak measurement theory a promising framework 
for the study of traversal times in quantum systems, a problem that does not fit well within 
standard measurement theory. 

In this paper, we show that the ambiguities which are present in the formalism when 
the traversal time problem is studied with the tools of standard measurement theory, |5j 
vanish in the framework of the weak measurement approach. However, the interpretation 
of the weak measurement results remains open. The outline of the paper is as follows: 
In Section 2 we present briefly the weak measurement theory (WMT), in a "minimalistic" 
way, i.e., concentrating on only those aspects of WMT that are directly relevant to the 
traversal time problem. We apply the technique to this problem in Sec. 3 and in Sec. 4 
show that several well known methods for defining and calculating average traversal times 
are particular realizations of the weak measurement approach. In Sec. 5 we go further and 
introduce an operator for the time spent in a region of space in an attempt to obtain higher 
moments of the traversal and dwell time distributions. A short discussion of open problems 
ends the paper. 

II. WEAK MEASUREMENT: A "MINIMALIST" FORMULATION 

In this section we describe the generic gedanken experiment and compare the standard 
measurement theory of von Neumann with the weak measurement theory of Aharonov and 
coworkers. For the scope of this paper we do not need to push the theory as far as Aharonov 
et al. and will limit the discussion to weak measurements on an ensemble of systems, 
staying clear of the more controversial issues of weak measurements on a single system and 
the reality of the wave function (i.e., the possibility of measuring the wave function of a 
single system). We use a minimalist approach to weak measurement theory treating it as 
a potentially useful extension of standard measurement theory, based on a "weak" system- 
apparatus interaction Hamiltonian. 

The experimental setup consists of a system E and a measuring device M evolving - 
when isolated - under the Hamiltonians and Hm, respectively. Let q be the canonical 
variable of the meter that we use as a pointer, and let 7r be its conjugate momentum. The 
corresponding operators are q and 7r with [tt, q] = —ih. 

To measure an observable A of the system E, let the system and apparatus interact 
through the Hamiltonian 

H mt = g(t)wA(t), (1) 
where g(t) = Gh(t), G is a constant and f_™ h(t)dt = 1. Let h(t) be non-zero only for 

te(ti,t f ). 

The system E and the meter M evolve independently with Hamiltonian Hq = Hy, + Hm 
until time t iy then undergo the interaction governed by H int , and, after time tf, continue 
their evolution under H . What is measured is the position of the meter at time tf. 

Let us denote by \ipo(t)), \cj) (t)), and | <^> (^) ) = iV'o(i)) ® \4>o{t)) the states representing 
the system E, the meter M, and their combination E plus M, respectively, evolving without 
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mutual interaction, and by the state of the combined system after the switching on of 

the interaction H int at time ij. Since the system E and the meter M do not interact before 
time t i: |$(t)) = \%(t)) for t < t { . 

For simplicity, we will consider Hm = 0, that is the state of the meter is static until the 
interaction is turned on, so that we can use = \4>o(ti)) for the state of the meter before 
time £j. Moreover, after the interaction is switched off, at tf, the state of the meter in each 
component of the linearly superposed entangled state no longer changes with time. 

In the Schrodinger picture, [0] 

\*(t f )) = &(t f ,t i )\*(t i )), (2) 
where U(tf,ti) is the evolution operator 

U(t f ,ti) = (ex P {-~£[H (t)+H int (t)]dty) , (3) 

and the +-subscript denotes time ordering of the integrals in the terms of the Taylor series 
expansion of the exponential function. In the following, we will indicate a state in the 
Heisenberg representation by omitting its dependence on time: for instance, |$) is the state 
!$(£)) in the Heisenberg representation, and is obtained as |$) = U(tf,t)$(t)). 



A. Standard Measurement 

In the von Neumann procedure, £/ tends to tj, i.e, h(t) ~ S(t — tf), and what is 
measured is the value of the observable A at the instant of time tf. 

In the time interval (ti,tf), Hi nt is the dominant term in the Hamiltonian and, from (0) 
and @, we have 

mt f )) ae-^/)^^)). (4) 
The probability density of pointer position q after the interaction is 

/(g) = (<%>(g|$> = £<$|an,g><an,?|S>, (5) 

n 

where {|a n (i))} is a complete set of eigenstates of A(t). Straightforward calculation yields 

f{q) = T,\ c n{t f )\ 2 \Uq-Ga n )\\ (6) 

n 

where c n {t) = (a n (t)\ip (t)) and <pi(q - Ga n ) = (q - Ga n \<pi). 

It is worth noticing that if the initial pointer position q is precisely defined, that is 
|0i(g)| 2 ~ 8(q), the probability density of the final position is a sum of quasi-delta functions 
in one-to-one correspondence with each of the eigenvalues of A. 
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1. Distribution of the pointer position 



The first two moments of the pointer position distribution are now easy to obtain. If we 
take an initial distribution of q centered at q — 0, the mean value of q at time tf is 

(q)f = mm = J g/(g)dg = G(A(t,)), (7) 

and the mean square value of q is 

(q 2 ) f = <<% 2 |$> = J q 2 f(q)dq = (q 2 ) t + G 2 (A 2 (t f )), (8) 

so that 

(A?/) 2 = (Q 2 )f - «<7>/) 2 = (^f + G 2 (AA f ) 2 , (9) 

where Aqf, Aqi, and AAf are the standard deviations of final and initial pointer positions, 
and of the observable A at time tf, respectively. The integrals without explicit limits are 
from — oo to +00. 



2. Verification of the unperturbed state 

It is interesting to calculate the probability that the state of the system S under obser- 
vation is not changed. In order to do so, we calculate the probability Pq of verification of 
the unperturbed state \ipo) at time tf, i.e. 



P Q (t f ) = ($|^o) (M®) = J (*\1> Q ,q)(il>Q,q\*)dq. (10) 
If we remember that |^o(*)) — J2n Cn{t)\a n (t)) we obtain 

Pa{tf) = E \cn(t f )\ 2 \c m (t f )\ 2 I <f>*{q - Gajfaiq - Ga m )dq, (11) 

n,m 

but if Aqi GAa, where Aa is the minimum difference between the eigenvalues of A 
(Aa = mm nj L m {\a n — a m \}), the integral in (pi]) is practically zero when n 7^ m, so that we 
can write 

Po(t f ) « E l^(t/)| 4 < max{\c n (t f )\ 2 }. (12) 



Equation (12) shows that the initial state is conserved only if it is an eigenstate of A; if 
this is not the case, the evolution of the system is strongly affected by the measurement. 
As will be shown in the next section, this problem does not exist in the weak measurement 
approach, due to the fact that the evolution of the system is perturbed only to order o(G) 
(by o(G) we mean a term such that lim G ^ o(G)/G = 0). 
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B. Weak measurement 



Weak measurement is characterized by the fact that the Hamiltonian for the interaction 
Hi n t is small enough to be considered as a small perturbation of the Hamiltonian Hq of the 
isolated system E, and the initial uncertainty in the position of the pointer q is much greater 
than G times the maximum separation between different eigenvalues of A. 

Most importantly, the interaction does not have to be impulsive, but can have a finite 
duration of time. This additional flexibility is a great advantage for measurements made 
over finite intervals of time. 

According to perturbation theory, || we can write 

\*(t f )) = |* (f/)> U (t f ,t)H int (t)\<S>(t))dt, (13) 

where 

Uo(t f ,t) = (e X p{-^iJ (t'K}) (14) 
is the evolution operator of the isolated system E. First order approximation on fll3"D gives 
|$(t,)) = (l + o(G))|$ (i/)) - 'j-Gtt J*' U (tf,t)A(t)\$ (t))h(t)dt, (15) 

where d(G) indicates a generic operator whose averages are o(G). 
If we introduce the hermitian operator in the Heisenberg picture 

X H (A) = f f U (tf,t)A(t)US(t f ,t)h(t)dt, (16) 

we can write (HD as 

|$) = [1 - j-GtiIh(A) + 6(G9po>. (17) 
Now we define A w , the weak value of the operator A, as 

A w = (%\i H (A)\<s> ) = f tf h(t)(Mt)\A(t)\Mt))dt (18) 

Jti 

The probability density of q after time tf is f(q) = and can be written, by using 

(0) and ([18]), as 

f(q) = ($ |1 + j-GttIh(A) + 6(G)|g>(g|l - ^-Gnl H (A) + o(G)|$ ), 

= ($ | exp(^-GA w 7r) + d(G)\q)(q\ exp(-^-GA w n) + 6(G) |$ ) 

= |^(g-G^)| 2 + o(G). (19) 

Except for terms of o(G), the final distribution of pointer positions is equal to the initial 
one translated by G times the weak value of A. It is worth noticing that if the interaction 
is impulsive (i.e., h(t) ~ S(t — tf)), we have A w w (A(tf)). 
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1. Distribution of pointer position 
The mean pointer position and the variance are, from ( |I7|) and (|19"D, respectively 

(q) f = = J qf(q)dq = GA W + o(G) (20) 

and 

(A?/) 2 = <<7 2 )/ " ((g)/) 2 = (^) 2 + o(G). (21) 

The average pointer position gives us the weak value of A; on the other hand, the variance 
does not give us additional information, because the weak measurement is very imprecise, 
due to the fact that the initial pointer distribution is very broad and the interaction is weak. 
Averaging over many identical experiments gives the right mean value, but does not tell us 
anything about the dispersion of the observed quantity, which is completely swamped by 
the dispersion in pointer position. 



2. Verification of the unperturbed state 

A fundamental property of a WM is that the evolution of E is practically not perturbed. 
In fact, verification of the initial state, using fllTD , yields 

P (t f ) = (®\i> )(M®) = l + o(G) (22) 

This means that several weak measurements of different observables on a single system can 
be performed. As a general property, and therefore even for non commuting observables, 
the order of successive measurements is not important. 



3. Conditional averages 

While conditional averages are not well defined within standard measurement theory 
||, they can be introduced in an unambiguous way within WMT, as a consequence of eq. 
(521) discussed above. Suppose that we want to measure the average of A conditioned to 
the verification of a given final state which is assumed, without loss of generality, to be a 
member \\ n ) of an orthonormal basis {|Xn)K for n = 1 . . .N, of the Hilbert space of E. Since 
|Xn)(Xn| and q commute, we can perform a standard measurement of both of them when the 
interaction is over, i.e., after time tf. Then, we keep only the readings of q corresponding 
to a positive verification of |xn)> an d calculate the "conditional" probability distribution of 
the collected readings f(qp n \ which is of the form 

f{q)(n) = < ^ix!><x^r > = Mq ~ G ^ n))|2 + (23) 



where 



s ^n\T H (Am = l_ [*' fr MA (t)\Mtm)dt. (24) 
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A^f> is the weak value of A for a system which is postselected in the state \x n ) (and 
preselected in the state |V>o))- To order o(G), the probability amplitude distribution of the 
meter's pointer is equal to the initial one translated by a quantity proportional to A$ . 
When defining A w in fllgp we did not specify a post-selected state; actually, to not perform 
a post-selection is equivalent to post-selecting the state |^o) because, as we showed in (p2|). 
verification of \ipo) is positive with probabililty very close to one. Therefore, Affl of (^) is 
equal to A w of ([18|) to order o(G) if |-0o) happens to be \Xn)- If it does not, then we can 
write \ipo) = J2nPn\Xn), where p n = (Xn\i>o), and have 

4 = EW (25) 

It is important to notice that, while t4„, is always real, Affl is in general complex valued. 

From (^|) and (^) we find that the conditional average and the standard deviation of 
the pointer position are, respectively: 

{q)f ] = / f(q) H qdq = GRe{A^} + o(G) (26) 

and 

(Aqff = (q^ - ((q)f^ = (Aq^ + o(G), (27) 

independent of n. 

In addition, from (^0]) and (|25|) we have 

(g)/ = El^| 2 (g)?' ) , (28) 

n 

that is, the well known sum law of conditional probabilities holds true for pointer position 
readings. 



III. WEAK MEASUREMENT AND TRAVERSAL TIMES 

Measurement of the time duration of some process requires that the observed system and 
the meter interact for a finite time, a situation for which the concept of weak measurement 
seems to be particularly well suited. Moreover, as we have just seen, WMT could also allow 
us to calculate conditional averages of a given temporal quantity for various outcomes of the 
unperturbed system. 

A well known and widely accepted result in the field of tunneling times is the dwell time, 
i.e. the average time spent by a particle in the region Q irrespective of its final state. If 
\ijjo) is the state describing the particle, the dwell time in the interval (tj,t/) is postulated 
to be 



(to) = r D (t h t f ) = f tf (Mt)\Pn\Mt))dt, (29) 



where Pq is the projection operator on the region Q. As can be seen, (p29]) is the mean 
value of Pq integrated over (tj,t/). It is hard to imagine this time as a result of a standard 



7 



measurement, because Pn is not a quantum non demolition (QND) variable |TT[] and, if 
h 7^ h, then Pn(ti) and Pnih) do not commute. 

However, (|29|) can be obtained as a result of a weak measurement. In fact, if we take 
A = Pq, and h(t) as constant in (ti,tf), the interaction Hamiltonian is 

H int = Gh(t)nP n , (30) 

and from (|TBD we have 



Pn w = 7 -/ (Mt)\Pn\Mt))dt. (31) 

t f - U 



Combining (p9|) and (|3~T|) yields the dwell time as 



<to> = MU,tf) = {t f - U)Pn w = hm (g)/( ^ (32) 

where we have used the fact that (g) / = GPq w + o{G). 

Suppose we are interested in the mean time spent in Q for some specified final state of 
the particle. Decomposition of dwell times in terms of particles evolving to a final state 
\Xn) is problematic within standard measurement theory, as has been pointed out many 
times: || the difficulty is that projection onto a region Q and projection onto a final state 
\Xn) involve non commuting operators, and there are no rules uniquely specifying how to 
build operators for quantities involving non commuting operators (this is also the reason for 
conditional probabilities being problematic). 

The ambiguity vanishes within the weak measurement approach: the weak value of Pq 
for a system postselected in the final state \x n ) is, according to 



= 1 I, w —. P ' (Xn(t f )\Pn\Mtf))dt (33) 

Therefore, the average time spent in Q from time tj to tf by a particle starting in the state 
|^o) and finally found in the state \\ n ) is 

(tn) (n) = r£> - fkzWdP = {tf - W PW}. (34) 

Summation over different final states holds: given |t/> ) = J2 n Pn\Xn) then, dropping the 
dependence on the time interval, we can write, from (|28|), (|32|), and 



(t D )=Y,\Pn\ 2 (tD) in) . (35) 



IV. WEAK MEASUREMENT AND WELL KNOWN METHODS FOR 
OBTAINING TRAVERSAL TIMES 

In this section we want to demonstrate that some well known approaches to the calcu- 
lation of tunneling times can be seen as particular examples of weak measurement, each 
corresponding to a different measuring apparatus. 
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In particular, we will focus our attention on methods based on the Larmor clock, |9| JT2|JT3 
on Feynman path-integrals, |15|,18|, and on absorption probabilities. [19] All of these proce- 
dures are based on the application of a small perturbation (a magnetic field, a real potential, 
an imaginary potential, respectively) to the region of interest. After that, the state of the 
particle evolves in time, and we attempt to extract the information about the time spent in 
the region of interest from some aspect of the perturbed wave function (the spin, the phase, 
or the amplitude, respectively depending on the kind of perturbation applied). In order not 
to perturb the evolution of the state too much, we let the perturbation tend to zero [18 . 
It has been demonstrated [^,^T|] that all the "probes" mentioned above lead to the same 
result. 

Let us now write two formulas that will be very useful in the remainder of this section. 
From Appendix A, the weak value of an operator A for a system postselected in the state 
\Xn), defined in ([24]), can be written as 



A (n) _ 9 (Xn(tf),7l\q\<5>) 
/i-,, — tttz ; — ; 



dG (Xn,vr|$ ) 



7T,G=0 



(36) 



Gtt=0 



where the second equality is true if q can be written as q = ihd/drr in the ^-representation 
and |$) depends only upon the product Gtc [as it obviously does for the interaction Hamil- 
tonian ([!])]. 



A. Real constant potential 



Let us start with a constant real potential applied only in Q and only for ti < t < tf. 
the perturbation Hamiltonian is Hi nt = Hy = f(t)VPn, with f(t) = 1 for t G (t{,tf) 



and zero otherwise. |18[ In order to translate this perturbation into the formalism of weak 
measurement, we can write V in the ^-representation as V = Gn/(tf — t{). Now the 
perturbative potential acting on the system S is of the form ( p0| ) . 

In this case, the weak value of the operator for a system postselected in the state 
\Xn) is, according to (|36|), 



/(Xn|r)(r,7r|^^|$)ci 3 r 



/(Xn|r)(r,7r|$ )d 3 r 



(37) 



Gtt=0 



We use the convention of omitting the limits of integration when the integrals run over 
the whole space. Given that V is proportional to tt, we can write $(r, V) = (r, 7r|<l>) and 
Xn(r) = (r\xn), so that (R3) becomes 



(to) 



(») 



(t f -U)Re{p£>} = Re 



(38) 



y=o- 



Note that (|58| ) is exactly the expression for the average time spent by a particle in the 
region Q obtained by using the Feynman path-integral technique. |15j If the final state is 
|r), i.e., the state corresponding to a particle found to be at r at time tf, the weak value of 
the average time is then 
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(tn) 



(r) 



Re 



ih d$(r,V) 
$(r, V) dV 



(39) 



v=o 



which is exactly the same expression obtained for the stay time defined in [18]. 



B. Pure imaginary potential 



A pure imaginary potential is often used in optics to simulate the absorption of photons 
by a material. What happens in this case is that the probability density of the particle is 
not conserved, because it decreases exponentially in Q, with a time constant proportional 
to the applied imaginary potential. The information about the average time spent in Q by 
the particle is therefore obtained by calculating how much of the total probability has been 
absorbed. 

The perturbation Hamiltonian in this case is PT 



H mt = H T = -f{t) l ^P a 



(40) 



which is of the form (|30|) if we put Y = 2iGit/(tf — £j). Analogously to (|37|) and (|38"D we 
have 



(to) 



(n) _ An) _ _ 



/x;(r)2ftj,$(r,r)d 3 r 



' D 



/Xn(r)$o(r, T)d?r 



(41) 



r=o 



where we have put $(r, Y) = (r,7r|$). This result, again, corresponds to the one obtained 
in [21]. 



C. Magnetic Field 



The well known Larmor clock method P~3J13|] involves applying an infinitesimal magnetic 



field in the z-direction, confined to the region fl. The spin, which is initially polarized in 
the x-direction, precesses in the x-y plane with the Larmor frequency ujl = eB/m when the 
spin is "in" Q. The spin polarization in the y-direction plays the role of pointer position. 
Let us consider as the perturbation Hamiltonian only the component which acts on the spin 
of the particle |2T 



H int = H B = f(t)^a z P n , (42) 

where a x , a y , and a z are the Pauli spin matrix operators. In this case 7r = ha z /2 acts as the 
pointer momentum and we put G = ojiitf — ti), so that fl42| ) takes the form (|30|). 

We have a x |^o) — l^o) because the initial state of the system is an eigenstate of a x . 
From 

h 

[&v, g^ll^o) = ift°x\i>o) = ift\^o) (43) 
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it immediately follows that q = a y and 7r = ha z /2 are the appropriate conjugate pointer 
operators. With this choice (^) becomes 



P, 



(n) 



Si a- 









8G (Xn,7T 





(44) 



7T,G=0 



and 



» 



r (n) 



Re 



9 (Xn,Tr\ 






8UJ L (Xn,7T 


l$o> 



(45) 



OJi,7T=0 



As is easy to see by comparison with (18) of Ref. [21], expression (f45|) for the time spent in 
VL is equal to the result obtained by Rybachenko IT2] and Baz'. p~3] 



V. HIGHER MOMENTS OF TIME DISTRIBUTIONS 



As is clear from (|2l|) weak measurements are not useful for obtaining higher moments of 
a distribution for the time spent in Q. In fact, the spread of final positions of the pointer is 
equal to the initial one to o(G). The only way within WMT of obtaining, say, the /th order 
moment of an operator A, is to build a meter sensitive to A 1 . This should have an interaction 
Hamiltonian of the form H^ t = Gh{t)^iA l (t) . In principle, there is no fundamental problem 
with this, and several meters can act simultaneously on the same system. 

The crucial point is that we need to use an operator for the time spent in Q, and not 
just the projector over Q as we did in section 3. In this section we will use the "sojourn 



time" operator previously introduced by Jaworski and Wardlaw. [22] It is consistent with 



the results of section 3 and 4, and is easy to obtain from the definition of mean dwell time 

A. An operator for the time spent in Q 

In the Heisenberg representation, the dwell time defined by ( p9| ) can be written as 

(to) = T D (ti,t f ) = (V>o|W|V>o> (46) 

if we just define 

t nH = f S U (t f ,f)P n U*(t f ,t')dt' = (t f -U)Z H (P a ), (47) 

J ti 

In the Schrodinger representation, the operator to corresponding to is 

in(t) = U*(t f , t)i nH Uo(tf, t). (48) 
For a gedanken experiment with a meter sensitive to tn, the interaction Hamiltonian is 

Hint = G x h{t)%M^ (49) 
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where h(t) = (tf — U)' 1 for t G (ti,tf), and otherwise; 7Ti and g\ are the conjugate 
momentum and position of the meter's pointer. From (48) if follows that Z#(£q) defined by 
(|j~6|) is equal to t^ H . Application of ( |I~8"D and ( |24] ) then leads to 



tQ W = {ipo\ZH(ia)\ipo) = (ipolinnlipo) , 



An) = (Xn\tnH\A) 
QW (Xn\tpo) 



(50) 



If we take (t D ) defined in (^), and (t^) defined in (|3~4"D, we can write 



(t D ) = lim 



7^— - \*d) - lim — p; — - ^e{t Qu , }. 

Lrl— >U Lr^ Gl^U (j-j 



(51) 



As can be seen, leads to the same result as Pa, in the measurement of average traversal 
times. 



B. Higher moments 

By the means of £q, we can measure any moment of order I of the distributions of times 
spent in Q. We need to use a meter whose corresponding interaction Hamiltonian is of the 
kind 

Hf\ t = GMttffa, (52) 

where tti and cji are the operators corresponding to the conjugate momentum and position 
of the meter's pointer. The average of the Zth power of the time spent in Q by a particle 
finally found in the state \x n ) is 

(n )^ 

(t l D )^= hm^ = Re{(4)W}, (53) 

with 

Only those pointer position readings corresponding to a postselected state \\ n ) are averaged. 
It is worth noticing that the sum rule of conditional averages is satisfied, i.e., if |?/>o) = 
E„Pn|Xn), then, for any integer I, 

(4) = Ebnl 2 (4) (n) (55) 

n 

It is also important to point out, while (t l D ) is positively defined, the conditional averages 
(t l D )^ are not. The lack of this important property has to prevent us from interpreting 
these quantities as the moments of a distribution of actual times spent by the electron in 
the region Q. 
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C. Comparison with some results in the literature 



The second moment of tjj, according to fl53|) and (p4|), is (t D ) = (^ol^jWI^o); if we 
remember that i^n = T H (pQ)(tf — ti), we obtain 

(tl} = (tf-U) 2 (HlUPnMo) 

= {tf-Uf J d 3 r(^ \l H (P n )\r)(r\I H (P n )\^ ) 

' d s r\t^\Mr,t f )\ 2 (56) 



where, as can be easily obtained from ( |3~3"D and (0), t$ is the weak value of the time spent 
in Q by a particle finally found in r. 

Eq. ([56]) is essentially equal to the result obtained for the second moment of the dwell 
time by a few works based on the path-integral approach. |]I6|- [T8|J25[ 

We would also point out that the second moment of the time spent in Q for a particle 
which is post-selected in position r at time tf, i.e., 

/f \(r) = Rc { ) _ (HPAh + *ng£#o) (5J) 

I Wo) J (MAl^o) 

where P r = |r)(r| is in general different from the time proposed in Ref. [18] on the 
basis of the path integral approach, that, in this formalism, would be equal to t 2 ^ = 
(^o \tn H PJn H | V»o) / (V>o I -Pr | ipo) 



D. Relation between higher moments and the measurement of the first moment 

In this section we show that the higher order moments of to obtained in Sec. 5.2 
can be obtained also from the wave function |<3>) of the system plus meter perturbed by 
the Hamiltonian for the first moment H int = Gih(t)ifit^i(t). In fact, if we multiply both 
numerator and denominator of fl54]) by (7r|0j), and substitute ( P8| ) in the numerator, we 
obtain 



(58) 



Gi=0 



If we put A = 7riGi, so that H mt = Xh(t)ia(t), and call $(A, r,t/) = (711, r|$(t/)) we can 
write for any integer / 

fe),r, = Re {i^)( ! 4)'* (A - r -«) < 59) 

Let us just point out that, while the form of (^) is exactly equal to the l-th complex 
moment of the dwell time distribution obtained on the basis of path integrals [|I6| , |TT|| , the 
meaning is substantially different, since the perturbative Hamiltonian used in path-integral 
approaches is of the kind H p \ = \h{t)P^i{tf — tj), while the perturbative Hamiltonian used 
for obtaining (p9|) is H^ t given by (fK]). It is clear, for example, that the former is local in 
space, while the latter is not. 
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VI. DISCUSSION 



Steinberg Q|| has argued that weak measurement theory is a promising tool for the 
study of the traversal time problem. Its major advantages over the standard measurement 
theory are the flexibility to treat interactions between a system and a measuring apparatus 
that are extended in time, and the possibilty of defining conditional averages for events 
corresponding to non commuting operators. Both these properties are due to the fact that 
a weak measurement prevents the wave function of the system from collapsing. 

We have shown that within WMT not only mean dwell and traversal times but also the 
averages of any higher powers of the time spent by particles in a region Q, conditioned to 
any final state of the system, can be mathematically defined in terms of the outcome of 
gedanken experiments. 

Unfortunately, there are severe problems of physical interpretation. As already pointed 
out for the special cases of the Larmor [I8,Z3| and Salecker-Wigner clocks |2?J, WMT may 
predict negative results for the average time spent by reflected particles on the far side of 
a barrier. In addition, as shown here, the conditional averages of any power of the time 
spent in fl are not positively defined within WMT. These unphysical results prevent us from 
interpreting them in terms of actual time spent by particles in the spatial region Q. 

To remain on firm ground, we are compelled to consider them as just quantities with 
the dimensions of time describing the response of a degree of freedom q of an apparatus to 
an interaction with particles that is constant in time over a finite time interval, linear, and 
proportional to a particle's presence in Q. Clearly, further investigation is required to learn 
whether these quantities can be fruitfully used to describe the time-dependent behaviour of 
E itself, i.e., apart from the particular interaction with the meter. 
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APPENDIX A: DERIVATION OF (36) 



We can start from the Eq. (p4[), where Affl is defined. If we multiply both numerator 
and denominator by (vr|0j) for n = we have 



(n) = (xn,7r|J H (i)|$ ) 
w <Xn,vr|$ ) 



(Al) 



7T = 



Now, we have just to remember that, 1 = [q, it]/ih and to substitute this formula into (|A1|) 
in order to obtain 
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/t( n ) 



1 

ifi 



(Xn,vr|g7TJH(A)|$ ) (x„, 7r|gJ H (A)|$ Q ) 

7T- 



(Xn,7r|^o) 



(Xn,7r|$o) 



7T=0 



the second term of this expression vanishes for it = 0. If we substitute 
the first term to the right of (|A2|) , we obtain Eq. (|36]). 



(A2) 

for Z = 1 into 



APPENDIX B: A FEW FORMULAS FROM PERTURBATION THEORY 



Let | and H^ t (t) be the system wave function and the interaction Hamiltonian, 
respectively, in the interaction representation, M i.e., 



\^(t)) = u (tf,t)mt)), 

Hg>(t) = Mtf,t)Hint(t)U;(t f ,t), 



(Bl) 
(B2) 



where Z7o(t/,t) = exp{—i/h f t f H (t')dt'} is the evolution operator. 

From ([Bip we have that |$/(t/)) = |$(t/)) = |$) and |$/(t<)) = |$o(*/)) = |$o), 



therefore 



(B3) 



where the +-subscript denotes time-ordering. 

If we take JEZj nt = Gh{t)7rA{t) as given by ([I]), with /i(t) = (t^ — tj)" 1 for t e (ti,tf) and 
zero otherwise, and put it in (fB2^) and (|B3), we obtain 



i*> = («p^^(i)n i$ ). 



Writing the exponential in 



m=0 



as a sum yields 
1 /G 



oo i , /~< \ m 

\*)= E-lb- vr m ([X H (i)r) J$ ), 



from which we obtain 



dG 7 



G=0 



y i (iMA)r) + \%), 



(ih) 



(B4) 



(B5) 



(B6) 



If we choose A(t) = in(t), we have the additional advantage that Th^q) = inn does not 
depend on time, so that time-ordering does not matter, and we can write 



dG l 



g= ( in >y 

from which we have, after projection onto the state \Xn,^), 



(B7) 



(B8) 



G=0 
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